Quantum calculation of vortices in the inner crust of neutron stars. 



O 
O 

(N 

X> 

PL, 
On 



> 
oo 

(N 
O 
(N 
O 
VO 
O 



o 

=5 



X 



P. Avogadro 0,6 , F. Barranco c , R.A. Broglia a,fc ^, E. Vigezzi'' 
a Dipartimento di Fisica, Universita degli Studi di Milano, via Celoria 16, 20133 Milano, Italy. 
b INFN, Sezione di Milano, via Celoria 16, 20133 Milano, Italy. 
c Departamento de Fisica Aplicada III, Escuela Superior de 
Ingenieros, Camino de los Descubrimientos s/n, 41092 Sevilla, Spain. 
The Niels Bohr Institute, University of Copenhagen, Blegdamsvcj 17, 2100 Copenhagen 0, Denmark. 

(Dated: February 9, 2008) 

We study, within a quantum mechanical framework based on self-consistent mean field theory, 
the interaction between a vortex and a nucleus immersed in a sea of free neutrons, a scenario 
representative of the inner crust of neutron stars. Quantal finite size effects force the vortex core 
outside the nucleus, influencing vortex pinning in an important way. 
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Neutron stars (pulsars) usually rotate with such a pre- 
cision that they are known as the best timekeepers in 
the universe. But every so often their rotation rate 
increases 0. Anderson and Itoh proposed that these 
glitches can be viewed as "vorticity jumps", equivalent 
to "flux jumps" in a superconducting magnet [2j. 

The density of the neutron star increases going from 
the surface to the interior, and when it becomes larger 
than about n = 5 xl0~ 4 fm ~ 3 , it becomes energetically 
favourable for some of the neutrons and all of the protons 
to lump together in extremely neutron rich Sn-like nuclei 
(that is nuclei containing 40-50 protons and about one 
hundred neutrons), surrounded by a sea of free neutrons. 
These nuclei form a pure electrostatic "Wigner" lattice. 
As one goes deeper into the inner crust, the lattice step 
decreases, going from about 90 fm at n — 5 x 10~ 4 fm -3 , 
to about 30 fm at n = 7x 10 -2 fm ~ 3 . There is strong evi- 
dence which testifies to the fact that in this density range, 
the neutrons are superfluid (cf. e.g. [3,4]). It has been ar- 
gued that the angular velocity of the superfluid in the in- 
ner crust of a neutron star changes either by vortex creep 
or by vorticity jumps, the latter causing the glitches. Al- 
though this scenario has been investigated in a number 
of publications [j| , there remain many open questions as- 
sociated with the variety of approximations used in these 
studies. In an effort to shed light into these questions, we 
report here the first fully quantal, self-consistent calcula- 
tion of a vortex, taking into account the inhomogeneous 
character typical of the inner crust of a neutron star. 

This allows one to calculate the pinning energy, namely 
the difference between the energy cost to create a vortex 
far from the nucleus and on top of it. Our calculations 
are performed for a single Wigner cell, considering a sin- 
gle nucleus, in keeping with previous estimates [f|, which 
showed that the distance between two neighbouring nu- 
clei of the Coulomb lattice is large compared with typical 
vortex dimensions, possibly except for the deepest layers 
of the inner crust. 

The calculations were performed solving the Hartree- 
Fock-Bogoliubov (HFB) equations (often called the De 



Gennes equations), within a cylindrical box of radius 30 
fm and height 40 fm, imposing that the wavefunctions 
vanish at the border of the cell. We have assumed that 
the vortex-nucleus system is axially symmetric ( with 
the vortex directed along the z— axis), but we have con- 
strained the proton density distribution (associated with 
deeply bound states) to have spherical symmetry. 
The De Gennes equations have the form 
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and have to be solved simultaneously with the number 
equation. 

The kinetic energy operator (including the effective 
mass associated with the Skyrme interaction) is denoted 
by K, V(p, z) being the self-consistent Hartree-Fock 
mean field and A the (S = 0) pairing field. It displays 
the functional form 



A(p,z,, 
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where p is the distance to this axis in the x — y plane and 
<t> is the azimuthal angle, while g denotes the density- 
dependent strength of a (contact) pairing interaction. 

Equations (l)-(2) allow for different solutions, which 
can be labeled by the vortex number v (y = 0,1,2...). 
We have mostly considered the v = 1 case, in which 
each Cooper pair carries one unit of angular momentum, 
its projection along the z— axis being +1. Note that for 
v = one recovers the Negele-Vautherin results |7J. 

The Skll force [g was used for determining the 
HF field, while in the pairing sector we have adopted 
the density-dependent contact interaction introduced in 
ref. 0, corresponding to the values g[n] — —481(1 — 
(n/n o ) 0A5 ) MeV fm 3 , where n a = 0.08 fixT 3 . This inter- 
action reproduces the values of the pairing gap calculated 
with the Gogny force in uniform neutron matter. 

The <fi dependence of the quasiparticle amplitudes U, V 
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is 

U a ( P , z, 0) = U a ( P , z)e a ^- V a ( P , z, 0) = V a (p, z)e^- v ^. 

. (3) 

The quasiparticle amplitudes are expanded on a basis of 
(free) single-particle wavefunctions inside the cylinder 

U a (p,z) = J2 U « m <l a (p)Xm(z), (4) 

nm 

and similarly for V a (replacing ip n ,i a with ip n ,i a -v)- 
The functions Xm(z) are (longitudinal) plane waves and 
V-Wc (p) (radial) Bessel functions, associated with a cylin- 
der with perfecly reflecting walls, l a being the single- 
particle angular momentum along the cylinder axis, cho- 
sen as quantization axis. To limit computational com- 
plexity we have not included the spin-orbit interaction 
term in the calculation of the single-particle levels. Con- 
cerning the protons, we have solved the v = equations 
in a spherical box of radius R p = 15 fm, using the same 
Skll force as for the neutrons. The effect of the neutrons 
on the (deeply bound) protons, has been included after 
a spherical average of the various neutron distributions 
has been carried out. 

In what follows we discuss the results associated with 
the value Ep = 5.8 MeV of the Fermi energy. The spatial 
dependence of the (neutron) density n and of the pair- 
ing gap A for a nucleus immersed in the neutron sea is 
shown in Fig. 1. The results for the density essentially 
coincide with the results of Negele and Vautherin for the 
corresponding (spherically symmetric) Wigner cell. The 
radius R of the nucleus is about 7.5 fm, the associated 
diffusivity being about 0.9 fm. Far from the nucleus the 
value of the pairing gap is about A U nif «2.2 MeV, equal 
to the value obtained in the case of uniform neutron mat- 
ter with the adopted pairing interaction, while at the nu- 
clear surface Awl MeV. The strong suppression of the 
pairing gap inside the nucleus is a consequence of the 
density dependence of the pairing interaction which re- 
flects the behaviour of the 1 S'o phase shift as a function 
of the relative kinetic energy of the pair of interacting 
nucleons. 

Let us now study the modifications induced on n and 
A by the presence of a vortex. We first discuss the case 
of free neutrons, a system which mimics, exception made 
at the edges of the cylinder, uniform neutron matter. In 
Fig. 3(a) we display the associated pairing gap for the 
case of a v — 1 vortex. It is seen that A vanishes along 
the z-axis. For small values of p, p < 3 — 4 fm, the gap 
increases linearly as a function of p. Defining the vortex 
core as the value of p for which A(p core ) = A un j//2, 
one obtains p core ~ 2 fm, a value which is similar to 
that of the correlation length £ m Tivf /2A(ps 5 fm). For 
larger values of p the gap increases more slowly, gradually 
approaching the value A un i f . 

Concerning the density distribution (see Fig. 2(a)), 
an axially symmetric depletion around the vortex axis is 
observed, with a radius (defined as the value of n at half 
saturation density) of the order of 3 fm, which is also of 



the order of the average distance between particles (r s ss 
1.92/fc F « 1.02/(£; F /20) 1 / 2 fm » 3.6 fm). The above 
results are similar to those found in ref. [10] within the 
framework of nuclear energy density functional, taking 
into account the different asymptotic values of the pairing 
gap (cf. also H3). 

We now turn our attention to the case of a v — 1 vortex 
pinned on the nucleus. The associated neutron density 
and pairing gap are shown in Fig. 2(b) and 3(b). Com- 
paring with the results displayed in Figs. 2(a) and 3(a) 
respectively, it is seen that both the density depletion 
around the vortex axis and the vortex core are strongly 
influenced by the presence of the nucleus. In fact, the 
vortex is seen to skate on the surface of the nucleus dis- 
playing a small penetration length. In other words, the 
vortex is essentially expelled by the nucleus. The pairing 
gap at the surface of the nucleus is strongly suppressed, 
both compared to the case of a isolated nucleus in the 
absence of the vortex (cf. Fig. 1(b)) and compared to 
the case of a vortex in uniform (neutron) matter (cf. Fig. 
3(a)). In particular, at p ~ 7 fm (and z =0) A displays 
a value of about 0.5 MeV (Fig. 3(b)), 1 MeV (Fig. 1(b)) 
and of 2 MeV (Fig. 3(a)), respectively. Furthermore, the 
presence of the nucleus delays the rise of the pairing gap 
to its asymptotic value by 5-7 fm as compared to the uni- 
form case (for example, in the z = plane, A becomes 
equal to 1.75 MeV at p ~ 11 fm, compared to p = 6 fm 
in the uniform case). 

There are two contrasting effects which influence the 
behaviour of a v = 1 vortex in the presence of a nucleus 
and lead to the results discussed above. First, the small 
(large) value of the pairing gap inside (outside) the nu- 
cleus favours pinning, because one saves pairing energy 
placing the vortex inside the nucleus. Second, to build 
a v = 1 vortex inside the nucleus requires the formation 
of Cooper pairs out of single-particle levels of opposite 
parity. This is strongly hindered by the spatial quantiza- 
tion associated with finite size nuclear effects which leads 
to a distribution of levels around the Fermi energy es- 
sentially displaying either positive or negative parity (in 
other words, a very small v = l(ir = (—1)" = —1) phase 
space [13 )• The consequences of spatial quantization on 
a v = 1 vortex can be further clarified by comparing the 
corrispondent solution of the De Gennes equations with 
that associated with a v = 2 vortex (see Figs. 3(c) and 
(d): note that within the present context the stability or 
less of this solution is immaterial). It is seen that the 
pairing gap in the case v = 2 displays a quadratic, rather 
than a linear dependence on p, reaching the asymptotic 
value at a larger distances from the z— axis than in the 
v = 1 case. It can also be seen that the pairing gap in 
the case v = 2 is modified very little by the presence of 
the nucleus. This is because v — 2 vortices thrieve on a 
subspace of single-particle levels all displaying essentially 
the same parity. 

We shall now consider the energy associated with the 
various situations shown above, for the case of a v = 1 
vortex, and calculate the correspondent pinning energy. 
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FIG. 1: Density and pairing gap calculated for a nucleus 
immersed in the neutron sea, at the Fermi energy Ef = 5.8 
MeV. They are spherically symmetric, but are calculated in 
the cylindrical box described in the text (radius 30 fm and 
height 40 fm). Only the region < z < 15 fm and < p < 15 
fm is shown. 
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FIG. 2: Density associated with av = l vortex, calculated 
in the cylindrical box described in the text (radius 30 fm and 
height 45 fm) , without and with a nucleus at the center of the 
box. Only the region < z < 15 fm, < p <15 fm is shown, 
(a) Density associated with a v — 1 vortex in the cell without 
the nucleus, (b) Density associated with a v — 1 vortex in 
the presence of the nucleus. 



In our calculation, the energy E to t of a given configura- 
tion receives contributions from three sources: the kinetic 
energy, the mean field potential (HF) energy, and the 
pairing energy, so that E tot = E kin + E pot + E pair . These 
contributions are displayed in Tabic 1 for the uniform 
case (no vortex, no nucleus in the cell: we call the value 
of the total energy Ejj), for an interstitial vortex (no nu- 
cleus in the cell, Ej), for an isolated nucleus (no vortex 
in the cell, En), and finally for a vortex pinned on the 
nucleus (Ep). For the first two cases, our calculation is 
essentially the same as the one performed in ref . ^(J , al- 
though using a different pairing interaction. The energy 
cost to create an interstitial vortex or a pinned vortex 
are given by Eju = Ej — Ejj and by Epjy = Ep — En, 
respectively In order to be meaningful, the cost to 
create a vortex should refer to two systems with the same 
number of neutrons moving in the same box. The cor- 
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FIG. 3: Pairing gap associated with a v = 1 and a v — 2 
vortex, calculated in the cylindrical box described in the text 
(radius 30 fm and height 40 fm), without and with a nucleus 
at the center of the box. Only the region < z < 15 fm, 
< p <15 fm is shown, (a) Gap associated with a v = 1 
vortex in the cell without the nucleus, (b) Gap associated 
with a v — 1 vortex in the presence of the nucleus, (c) Gap 
associated with a v — 2 vortex in the cell without the nucleus, 
(d) Gap associated with a v — 2 vortex in the presence the 
nucleus. 







Ekin 


Epot 


Epair 


Etot 






Uniform 


6841.9 


-1735.3 


-1322.1 


3784.5 






Vortex, int. 


6776.0 


-1737.5 


-1203.9 


3834.6 






Nucleus 


9971.9 


-5784.0 


-1274.5 


2913.4 






Vortex, pinn. 


9893.4 


-5806.1 


-1120.5 


2966.8 





TABLE I: The total energy E to t, subdivided into the three 
contributions arising from the kinetic energy Ekin, the poten- 
tial energy E pot and the pairing energy E pa i r , is shown for 
each of the four configurations discussed in the text (Wigner 
cell without the nucleus, cell without the nucleus and with an 
interstitial vortex, cell with the nucleus, cell with the nucleus 
and a pinned vortex). 

responding solution of the self-consistent equations (1-4) 
leads to a slight change 6Ep of the Fermi energy in go- 
ing from a system without a vortex to a system with a 
vortex. For an infinitely large box, the change tends to 
zero. In the present calculation the change was SEp w 
0.05 MeV. 

We have evaluated the error associated with the fi- 
nite mesh size used in integrating these equations (Ap = 
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Az = 0.25 fm). In particular we have found that the 
value of the pinning energy, although obtained from the 
subtraction of large numbers, is remarkably stable with 
respect to changes in the box size. In fact we estimate 
that the absolute value of the error associated with this 
quantity is less than 2 MeV. 

The main contribution to Eiu = Ei — Ejj as reported 
in Table 1 originates from the pairing energy (AE pa i r = 
118.2 MeV), and is associated with the decrease of the 
pairing field in the region close to the vortex axis (see 
Fig. 3(a)). At the same time, there is a reduction of 
kinetic energy due to the reduced population of levels 
above the Fermi energy. On the other hand, one also ex- 
pects a positive contribution to the kinetic energy in the 
presence of a vortex, associated with the irrotational flow 
around the axis. In the present case the balance turns 
out to be negative (AEkin = —65.9 MeV). The third con- 
tribution to Ejxj arises from a modification in the (HF) 
mean field. This energy is related to the redistribution of 
the particles around the vortex, which, in the case under 
discussion, is small (AE pot ~ -2.2 MeV). Adding the three 
contributions leads to Em = AE pa i r + AEkin + AE pot = 
118.2 - 65.9 - 2.2 MeV = 50.1 MeV. Analogous consider- 
ations can be made to calculate the energy cost Epn to 
create a vortex pinned on the nucleus. We have in this 
case E PN = AE pair + AE kin + AE mean = 154.0 - 78.5 
-22.1 MeV = 53.4 MeV. It is remarkable that the cost 



in pairing energy for a vortex pinned on a nucleus (154.0 
MeV) is larger than the cost involved in the creation of a 
vortex in uniform matter (118.2 MeV). One might have 
thought that because the vortex tends to lower the pair- 
ing gap it would cost less pairing energy to create it in 
the region occupied by the nucleus, where the pairing gap 
is lower than in uniform matter (that is, far away from 
the nucleus, see Fig. 1(b)). However, because the vortex 
tends to avoid the nucleus, its presence affects the neu- 
tron pairing gap essentially only at the nuclear surface. 
For example, for z w 0, p w 7 fm in the uniform case, the 
pairing gap changes from A un if (2.2 MeV, v — situa- 
tion ) to a value of 1.7 MeV ( v = 1 vortex, cf. Fig. 3(a)). 
In the pinned case the pairing gap changes from 1.6 MeV 
(v = situation) to about MeV (v = 1 vortex). In 
other words, the main effect of a vortex in the pinned 
situation is to zeroth the pairing gap at the nuclear sur- 
face, an effect involving an important amount of pairing 
energy. The pinning energy is defined as the difference 
between the energy cost to create a pinned vortex, and 
to create a vortex in uniform matter @. For the Fermi 
energy under consideration (Ep = 5.8 MeV) it amounts 
to E PN - E w = 3.3 MeV [13]. 

One can conclude that quantal finite size effects, in par- 
ticular the spatial quantization leading to shell structure, 
have important effects on the vortex-nucleus interplay. 
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